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ABSTRACT 

We report a discovery of proto-cluster candidates around high redshift radio galaxies 
at z ~ 2.5 on the basis of clear statistical excess of colour-selected galaxies around 
C***~ , them seen in the deep near-infrared imaging data obtained with CISCO on Subaru 

Telescope. We have observed six targets, all at similar redshifts at z ~ 2.5, and our 
data reach to J = 23.5, H = 22.6 and K = 21.8 (5a) and cover a 1.6' x 1.6' field 
centered on each radio galaxy. We apply colour cuts in JHKs in order to exclusively 
search for galaxies located at high redshifts, z > 2. Over the magnitude range of 
19.5 < K < 21.5 we see a significant excess of red galaxies with J — K > 2.3 by a 
factor of two around the combined radio galaxies fields compared to those found in the 
i-C ■ general field of GOODS South. The excess of galaxies around the radio galaxies fields 

becomes more than factor of three around 19.5 < K < 20.5 when the two-colours cuts 
q ■ are applied with JHKs. Such overdensity of the colour-selected galaxies suggest that 

!— i ' those fields tend to host high density regions at high redshifts, although there seems to 

\ be the variety of the density of the colour-selected galaxies in each field. In particular, 

■ two radio galaxies fields out of the six observed fields show very strong density excess 

and these are likely to be proto-clusters associated to the radio galaxies which would 
evolve into rich clusters of galaxies dominated by old passively evolving galaxies. 

5_j ■ Key words: galaxies: clusters: — galaxies: evolution: — galaxies: high-redshift 



1 INTRODUCTION 

From many previous studies based on fundamental planes 
and colour-magnitude (CM) relations in clusters up to 

2 ~ 1.3, it appears that major star formation in the 
brightest cl uster galaxies was virtual l y completed by 
z ~ 2 (e.g.. iBower. Lucev. fc Ellis! Il992t lEllis et al.lll997t 

fegdam^^t^JJ^^^^^^^^^2I^^^^^^^^^^SI^^^g3 

llQQSt Ivan Dokkum et al.l ll99Si iBlakeslee et alJ I2003TT 

Therefore, one expects to find proto-clusters dominated 
by already passively evolving galaxies even beyond z ~ 2. 
Searching for such ancestors of present-day rich clusters 
of galaxies is not an easy task, however, since it requires 
a large volume to survey in order to find such rare ob- 
jects. Radio galaxies can be viewed as good markers of 
such high density regions at high redshifts, since they 
are among the most massive galaxies at any redshift 
(M st ar > 10 11 M m . lR"occa-Volmerange et alJl2004T) . In fact, 
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about half of powerful radio galaxi es at intermediate 
redshift inhabit rich environments (Hill fc Lillvl Il99lt 
iBornancini. Lambas. fc De Breuckl 120061) . The hypothesis 
that radio galaxies point us towards proto-clusters at high 
redshift has been further supported by the recent finding 
of a correlation betw een AGN activ i ty an d bulge mass 
(Magorrian relation, iGebhardt et alJ l2000f) . Given this 
relation, it is natural to imagine that massive cD galaxies 
sitting in deep potential well of clusters tend to host massive 
powerful AGN and hence they are viewed as radio galaxies 
at high redshift. 

Narrow-band surveys of Lya emitters around high-z ra- 
dio galaxies have been intensively conducted by the L eiden 
group (e.g.. iKurk et al.ll200l IVenemans et alJ 12002^1 over 
many years who discovered a large number of emitters at 
the same redshifts of the radio galaxies, and a large frac- 
tion of those emitters have been spectroscopica lly confirmed 
(e.g.. lKurk et alJl200l [Venemans et al.ll2002l) . These over- 
dense regions of star forming galaxies are embedded in large 
scale structures at high redshifts and are likely to evolve 
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Table 1. Summary of the observational data 



Name 


Redshift 


J 


H 


Kp 


PSF 






(sec) 


(sec) 


(sec) 


(arcsec) 


4C -00.62 


2.527 


6560 


2220 


3880 


0.67 


4C -00.54 


2.363 


5760 


3200 


3200 


0.65 


4C +23.56 


2.483 


6800 


4525 


2640 


0.67 


MRC 2104-242 


2.419 


5400 


3630 


2880 


0.60 


4C +04.81 


2.594 


5400 


4320 


2880 


0.45 


MRC 0406-244 


2.440 


5400 


3960 


2560 


0.45 



into massive systems at the present day, such as clusters of 
galaxies. 

Alternative approach is to find evolved populations at 
longer wavelength. This i s complementary to the emitter 
surveys. iBest et all <l2003l) found an excess of red galax- 
ies with R — K > 4 around radio-loud AGNs at z ~ 1.6. 
lHall et al] il998ft and iKodama fc Bowed ]2003l) also found 
statistical overdensities of red galaxies aroun d radio-loud 
AGNs at 1 g z g 1.5 in near-infrared surveys. IWold et al] 
(2003) also report overdensities of red galaxies in optical- 
NIR colours around 13 QSO's at 2 g z & 3. These studies 
indicate that there is also an excess of evolved population 
around the high-z radio galaxies, as expected if they are 
embedded in proto-clusters. 

We take the second approach and extend the previous 
analyses to higher redshift based on the near-infrared (NIR) 
imaging. And in this Letter, we report a discovery of the 
overdensity of NIR-selected galaxies in the fields of six radio 
galaxies at z ~ 2.5. 

The Vega-referred magnitude system is used throughout 
the paper. 



2 OBSERVATION AND DATA REDUCTION 

We performed the J// /("-bands imaging observations of the 
fields of the five radio galaxies an d one radio-loud QSO (4C 
+04.81) atz~ 2.5 with CISCO jMotohara et aljl2002D on 
the Subaru Telescope on 2005 March 21-22 and Septem- 
ber 19-20. CISCO has the field of view of about 1.8 x 1.8 
arcmin 2 with 0.105" pixel scale. The target information and 
exposure time are given in Tabled The weather conditions 
were stable during the observations, and seeing sizes were 
between 0.3" and 0.7" (FWHM). 

The data were r educed using the IRAF soft ware pack- 
age as described in iKaiisawa fc Yamadal J2005) . We per- 
forme d flat-fielding with the superflat frames for the CISCO 
data jMotohara et al]l2002t) . and subtracted the dark and 
sky background. Then the data were co-registered and com- 
bined. We also reduced the archival CISCO data for 4C - 
00.62 (5600 sec in J-band and 1040 sec in if-band) and 4C 
+23.56 (5600 sec in J-band, 1600 sec in //-band and 1200 
sec in i^-band) and combined together. These are included 
in Table Our final multi-band combined images have the 
slightly smaller field of view of about 1.6 x 1.6 arcmin 2 , 
which corresponds to ~ 0.77/i7oMpc on a side at z ~ 2.5 
(We assume a cosmology of Qo = 0.3, Qa = 0.7). For each 
target, the images were convolved with a Gaussian kernel to 
match the PSFs to the worst one. The final PSF for each 
target is also given in Tabled 



Source detection was performed in the K-band 
images using the SEx tractor image analysis package 
jBertin fc Arnoutsl ll996V We adopted MAG.AUTO from 
the SExtractor as the total /("-band magnitude of de- 
tected objects. For the colour measurements, we used the 
MAG_AUTO with the Kron factor of k = 2, which corre- 
sponds to the 20% smaller aperture size than that used in 
total magnitude with default value k = 2.5. Kron radius of 
each object was determined in the AT-band images and the 
same radii were used for all J, H, and /f-bands. 

The UKIRT Faint Standards were used for the flux cal- 
ibration. For the data obtained in the March run, we had 
to shift the //-band zero-point by 0. 13 mag in order to re- 
prod uce the stellar tracks in JHK llKideer fc Martm-Luisl 
2003) using the observed stellar objects in our combined im- 
ages. 

The 5a limiting magnitude at 1.2 arcsec diameter is 
about J = 23.5, H = 22.6 and K = 21.8, respec- 
tively. All the magnitudes are corrected for the Galactic 
extinction, which is estimated at the positions of the ra- 
dio igah5jd^s_(near_tlw_ce^^ field) based 
on lSchleeel. Finkbeiner. fc Davis] <ll998l) . 



3 COLOUR SELECTION OF HIGH-Z 
GALAXIES 

Figure shows J — K vs K colour- magnitude diagrams of 
the six observed radio galaxies fields. Type of symbols differ- 
enciate the galaxies selected with different colour cuts on the 
J — K vs. H — K diagram (Figure^ see below). The square 
in each diagram indicates the targeted radio galaxy (radio- 
loud QSO) at 2 ~ 2.5. Some of the fields contain a bunch of 
very red galaxies (J — K 2.3), in particular 4C —00.62 and 
4C +23.56 show sequences of red galaxies which are likely to 
be progenitors of the colour-magnitude relation of early-type 
galaxies commonly seen in lower redshift clusters at z g 1. 
Therefore these two fields are most plausible high- 2 cluster 
candidates. For reference, we also plotted an expected C-M 
relation at z — 2.5 using the GALAXEV popu lation synthe- 
sis spectral library feruzual fc Charlotll2003T) . We assumed 
a single burst with a duration of 0.1 Gyr and the formation 
redshift of 2f rm = 5. The observed red galaxies tend to be 
relatively bluer than the model prediction. In contrast, the 
MRC 0406-244 field contain few red galaxies, which may in- 
dicate the variety of the number densities of red companions 
around radio-loud AGNs at 2 ~ 2.5. 

In order to efficiently pick up plausible proto-cluster 
members associated to the central radio galaxies at 2 ~ 2.5, 
we make use of our multi-colour photometric data at NIR. 
We apply three different colour cuts as shown in Figure a 
combined colour-colour diagram in J — K vs. H — K for all 
the six fields. 

First is a simple and single colour cut of 

J-K>2.3 (1) 

as shown in the horizontal dotted line. This colour-cut is 
the same as is used for the selection of the Distant Red 
Galaxies at 2 £ 2, and many of thus selected galaxies are 
spectroscopically confirmed. (e.g.,|van Dok kum et aljl200l 
iForster Schreiber et al] I2004L iReddv et al]l2005l) . In fact, 
such red colours of galaxies can only be reproduced either by 
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Figure 1. Colour-magnitude diagrams of the six observed radio galaxies fields. Type of the symbols diffcrcneiate the adopted colour 
selections as indicated in Figure 2 (see also text). Dashed lines show the 2er limit at J-band. For objects with J-band fluxes lower 
than the limit, the lower limits of J — K color are plotted by symbols with arrows. A square in each panel indicates the targeted radio 
galaxy. Solid lines represent the expected C-M relation at z = 2.5 calculated from a GALAXEV single burst model with a duration 
of O.lGyr with Zf nrrn = 5. T he metallicity of the model is adjusted so that the model coincide with the C-M relation of Coma cluster 
jBower. Lucev. fc Ellislll992D at z = 0.02. 



Balmer/4000A-break galaxies with old populations at z £ 2 
or the foreground galaxies with heavy dust extinction. 
Second colour cut uses two colours as: 

J-K > 2 x (H- K) + 0.5 & J- K > 1.5 (2) 

which are shown by thick solid lines in the figure. Here 
we overplot the model colour tracks with various star for- 
mation histories using the GALAXEV library. As shown, 
the top left corner of the diagram separated by Equation 
(2) should be exclusively populated by galaxies located at 



2 <; z % 3. This criterion is met while the Balmer/4000A- 
break of galaxies falls between J-band and i/-band, while 
most of the foreground and background galaxies would not 
meet such criterion. Therefore, by applying such colour-cut 
in JHK, we can effectively achieve strong contrast of plausi- 
ble proto-cluster members at z ~ 2.5 against numerous fore- 
ground/background contamination. This two-colour-based 
selection also has a significant advantage over the first one 
based on the single J — K colour, since the two-colours 
cut can include younger galaxies or star forming galaxies 
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Figure 2. Combined J — K vs. H — K colour-colour diagram of all the six radio galaxies fields. The horizontal dotted line and the 
thick solid lines show the boundaries of our single colour selection of J — K > 2.3 and the two-colours selection with JHK, respectively. 
Different symbols (solid circle, open circle, triangle, cross) correspond to our different colour selections as shown in the legend (see also 
text). The size of the symbols for the colour-selected galaxies is scaled according to apparent magnitudes in _R"-band (the bigger, the 
brighter). The six open squares indicate the targeted radio galaxies at z ~ 2.5. The long arrow at the left side of the figure shows a 
reddening vector of Calzetti et al.(2000)'s law of E(£? — V) = 0.3 at z = 2.5. Dashed, dotted and dotted-dash curves represent the 
GAL AXE V models of evolutionary sequences over < z < 4 of galaxies with Zf OIm = 5 for various star formation histories. Pentagons 
show the points of z = 2.5 in the model sequences and shaded squares show those of z =1 ,2,3, and 4. 



at z ~ 2.5 which have relatively blue colours in J — K and 
would have been missed by the single J — K > 2.3 cut. It 
should be also noted that this selection is robust against 
dust extinction, since the reddening vector (the long arrow 
in Figure 1^1 is almost parallel to the bounda ry line of the 
J i/jf -selection. We also confirmed with the iKnapp et alJ 
(2004)'s star catalogue that most of cool M, L and T dwarfs 
should be excluded by this criterion, since they are either too 
blue in J — K or located on the right side of the boundary 
line of the JH Tf-selection. 

Different types of the symbols correspond to different 



colour selections. Filled circles show the galaxies which sat- 
isfy the both criteria simultaneously. Open circles represent 
the objects which meet only the second criterion, and open 
triangle indicate those which meet only the first one (i.e., 
DRGs which do not satisfy the second colour-cut hence are 
likely to be background red objects). All the others are 
shown by crosses. 

It is notable that all the five radio galaxies at z ~ 2.5 
(squares in Figure |2J satisfy the two-colours cut (the re- 
maining square outside of the top left corner is the radio- 
loud QSO 4C +04.81). This supports that our colour selec- 
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Table 2. Number densities of the colour-selected galaxies with 
K < 21.8 (arcmin -2 ) 



Field 


J-X>2.3 


.////•("-selected 


J///i"-selected 
& J-X>2.3 


4C -00.62 
4C -00.54 
4C +23.56 
MRC 2104-242 
4C +04.81 
MRC 0406-244 


2.28 ± 0.93 
2.78 ± 1.05 
5.53 ± 1.48 
1.15 ± 0.66 
2.31 ± 0.95 
0.38 ± 0.38 


6.08 ± 1.52 
2.38 ± 0.97 
6.72 ± 1.63 
2.30 ± 0.94 
1.16 ± 0.67 
0.38 ± 0.38 


1.52 ± 0.76 
1.19 ± 0.69 
3.95 ± 1.25 
0.38 ± 0.38 
1.16 ± 0.67 
0.00 ± 0.00 


-0062 & +23.56 
the other 4 fields 


3.87 ± 0.87 
1.68 ± 0.39 


6.97 ± 1.16 
1.68 ± 0.39 


2.71 ± 0.72 
0.65 ± 0.25 


all 6 fields 


2.39 ± 0.39 


3.39 ± 0.46 


1.32 ± 0.29 


GOODS-S 


1.66 ± 0.97 


1.75 ± 1.03 


0.66 ± 0.60 



tion works reasonably well, even though the strong emission 
lines from the AGN component could affect their colours 
fllwamuro et alJl2003h . 

Tabled list the number densities of the colour-selected 
galaxies with K < 21.8 in each radio galaxy field. Poisson 
errors of these densities are also shown. In addition, we used 
the public GOODS-South data (see next section for the de- 
tail) to estimate the degree of the field-to-field variance for 
the observed area of each radio galaxy field. A random field 
of 1.6x1.6 arcmin 2 was selected from GOODS-S field and 
the numbers of color-selected galaxies were counted. We re- 
peated this procedure and measured the one sigma variance 
in the numbers. Such field-to-field variance is shown as er- 
rors in the number densities for the GOODS-S field in Table 
121 and Figure |3] (see next section). 

It is clear that the two fields, 4C -00.62 and 4C +23.56, 
contain quite high number densities of galaxies that sat- 
isfy the two-colours cut (Eq. 2) and they are among the 
most plausible proto-cluster candidates. In particular, the 
4C +2 3.56 field has numerous red g alaxies with J — K> 2.3. 
In fact. lKnopp fc Chambers! (Il997tl also reported a high den- 
sity of red galaxies with I—K > 4 in the central 1.25 arcmin 2 
field around 4C +23.56. On the other hand, MRC 0406-244 
field has few colour-selected galaxy. There seems to be the 
variety of the density of the colour-selected galaxies in each 
radio galaxy field. In the following section, we compare those 
number densities of these colour-selected galaxies in the ra- 
dio galaxies fields with those in a general blank field. 



4 OVERDENSITIES OF HIGH-Z GALAXIES 

In order to quantify the overdensities of the colour-selected 
galaxies around the hig h-z radio gal axies, we use GOODS- 
South survey IIGiavalisco et al.ll2004fl as a comparison field 
for which the public VLT/ISAAC verl.5 data are available. 

The JHK-b&nds data of the GOODS-S has similar 
depths and similar PSF to those of our CISCO data, and 
are suitable for our purpose. The area used in the analysis 
is 94 arcmin 2 , which is significantly larger than that of our 
radio galaxies fields. We performed the source detection and 
photometry exactly in the same manner as for the CISCO 
images. 

In Figured! we compare the cumulative K-b&nd number 



densities of the colour-selected galaxies in the radio galaxies 
fields and in the GOODS-S field. In the left panel, we plot 
the number counts of all the Jf-selected objects and the 
galaxies with J — K > 2.3. Note that the excess of all the K- 
selected objects at the bright end (K <J 19) is due to stars in 
the radio galaxies fields, especially in the case of 4C +23.56 
field which contain numerous stellar objects (Figure^. The 
number counts of galaxies that are not selected by any color 
cut ( "the others" in Figure^ in the radio galaxies fields and 
GOODS-S are consistent, if we exclude these star excess in 
the radio galaxies fields by the selection such as J — K > 1.0 
or K > 19. The excess of all the /("-selected objects at K > 
19 is caused by the overdensity of the color-selected galaxies 
as seen below. 

The cumulative number density of the galaxies with 
J — K > 2.3 in all the six radio galaxies fields significantly 
exceeds that in the GOODS-S at K > 19.5. The former den- 
sity is twice higher than the latter over the magnitude range 
of 19.5 < K < 21.5. 

The middle panel of Figure [3] shows the cases for galax- 
ies selected by the two-colours cut in JHK (open circles). 
The galaxies shown by the filled circles have been further 
selected as having red colours of J — K > 2.3 in addition to 
the two-colours cut. We see a significant excess in densities 
of the J//A"-selected galaxies around the radio galaxies. In 
fact, the overdensity is more than factor of three compared 
to that of the GOODS-S at 19.5 < K < 20.5. Red subsam- 
ple with J — K> 2.3 also shows similar overdensity of about 
factor three over the range of 19.5 < K < 22. 

In the right panel of Figure |3] we show the same 
plots as in the middle panel but for only the two highly 
plausible proto-cluster candidates, namely 4C —00.62 and 
4C +23.56 fields. The overdensity of galaxies selected by 
the two-colours cut is now more than factor of five com- 
pared to the general field. The excess is even higher for the 
red subsample with J — K > 2.3. Such strong excess of 
NIR-selected (approximately stellar mass-selected) galaxies 
indicates that these fields are likely to be proto-clusters with 
evolved populations associated to the central radio galaxies 
at z ~ 2.5, which would evolve into rich clusters of galaxies 
today dominated by old passively evolving galaxies. In fact, 
these two fields mainly cause the significant excess seen in 
all the six fields, while the average density of the other four 
fields is consistent with that of the general field within the 
uncertainty (see Table 

In Figure |1J we show the spatial distribution of the 
colour-selected galaxies in the two proto-cluster candidates 
around the z ~ 2.5 radio galaxies. Despite of poor statis- 
tics, it is suggestive that those colour-selected member can- 
didates are distributed non-uniformly around the central ra- 
dio galaxies. Such non-uniform distribution of gala xies in 
proto- clusters has also been reported recently bv lCroft et alJ 
(2005) who found a filamentary structure in the spatial dis- 
tribution of spectroscopically confirmed Lya emitters asso- 
ciated to a radio galaxy at z = 2.16. The mass-selected 
members of proto-clusters may similarly trace such filamen- 
tary structures at least over the scale of ~ 0.8/i7oMpc at this 
observed epoch. 




Figure 3. Combined cumulative number counts of galaxies in the radio galaxies fields (solid lines) and those in a controle field GOODS- 
South (dashed lines). Left: For all the i^-selected objects (squares) and for red galaxies with J — K > 2.3 (triangles). Middle: For 
galaxies selected by two-colours cut (open circles) and for galaxies which satisfy the two criteria simultaneously (solid circles). Right: 
Same as the middle panel but only for the two most probable proto-clusters of 4C —00.62 and 4C +23.56. Errorbars on the radio galaxies 
fields are Poissonian (i.e., square roots of the observed numbers), and errorbars on the GOODS-S data indicate the field-to-field variance 
scaled to the same area as the radio galaxies fields (six, six and two times 1.6 X 1.6 arcmin 2 for the left, middle and the right panel, 
respectively). 
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Figure 4. Spatial distribution of the two most probable proto-clusters, 4C -00.62 and 4C +23.56. Symbols are the same as Figure l2l 



5 SUMMARY 

We searched for proto-clusters around the radio galaxies 
at z ~ 2.5 based on the deep JHK -bands imaging taken 
with Subaru telescope. In order to largely subtract fore- 
ground/background contamination and have high contrast 
of proto-cluster member candidates associated to the radio 
galaxies, we applied several colour selections; (1) a widely 
used single colour cut of J — K > 2.3 and (2) our new two- 
colours cut with JHK-haxiAs. In both cases, we see a clear 
excess in the number densities of colour-selected galaxies in 
the observed radio galaxies fields compared to the general 
field GOODS-South. The overdensity in the six combined 
radio galaxies fields is about factor of two for the galaxies 
with J — K > 2.3 and is factor of three for those selected by 
the two-colours cut, although the variety of the overdensity 



in each field is also seen. In particular, two radio galaxies 
fields, namely 4C -00.62 and 4C +23.56, show the highest 
overdensities of the colour-selected galaxies of more than 
factor of five. Therefore, we are likely to be witnessing the 
assembly and formation of proto-clusters with many evolved 
galaxies at z ~ 2.5 associated to the radio galaxies. 
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